Morphology— transitions at vicinal Cu— surfaces based on entropic step— step interaction 

and diffusion along steps 



Heike Emmerich 

Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, D 01187 Dresden, Germany 



o 
o 



An extension of the Burton-Cabrera-Frank model [Phil. Trans. R. Soc. London, Ser. A 
243, 299 (1951)] including diffusion along steps and entropic step-step interaction is introduced. 
This extended model is successfully applied to simulate experiments done in the group of H.-J. 
Ernst at CEA Saclay, France, for vicinal Cu surfaces. In particular, the rise of two qualitatively 
different morphologies can be explained by the competition of growth directions implied by the four 
distinct driving and restoring forces of the model. Each of those forces sets a length scale for the 
dominant mode of a developing step instability. This explains that the wavelengthes of meandering- 
instabilities observed at vicinal surfaces experimentally are usually larger than the one predicted by 
Bales and Zangwill [Phys. Rev. B 41, 5500 (1990)] theoretically. 
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During molecular beam epitaxy (MBE) appropriate 
conditions for the controlled growth of vicinal surfaces 
can be fine-tuned. This way one tries to fabricate either 
atomistically flat or nanostructured surfaces. The ability 
to exert control on structuring along the growth direction 
and fabricate substrates with smallest-scale-built-in pe- 
riodicies normal to the surface is well advanced. Efforts 
to still enhance abilities to determine functional proper- 
ties of a grown substrate now focus on the lateral struc- 
turing within one layer of growth. One direction along 
this line is to make use of the inherent instabilities due 
to the dynamics of the growth process itself Q. A nec- 
essary first step is to understand the basic wavelengthes 
of those inherent instabilities just as well as the kind of 
morphologies which will develop. 

More then ten years ago, Bales and Zangwill 1^ pre- 
dicted that a growing vicinal surface should undergo a 
step meandering instability when kinetic step edge barri- 
ers suppress the attachments of atoms to decending steps 
1^. According to their analysis, a straight step is linearly 
unstable against perturbations with wavelengths larger 
then Ac — 2Tr(2TL^) and a fastest growing wavelength 
at A„ = V2Xc (r = La = -^). Here F de- 

notes the deposition rate, D the diffusion constant for 
diffusion on the terraces and Cf,q the equilibrium concen- 
tration, 7 the step stiffness, fl the atomic area and k^T 
the termal energy. Even though the meandering insta- 
bility has meanwhile been observed in experiments and 
simulations Q , the quantitative prediction of the Bales- 
Zangwill analysis could not be recovered in many of those 
experimental findings. This point received much interest 
in view of the recent experiments by Maroutian et al. 

at CEA Saclay. Analytical efforts to resolve the dis- 
agreement between experimental measurements and the- 
oretical prediction led to a precise study of the limit of 
desorptionless growth. For to derive a single evolution 
equation in the weakly nonlinear regime this case has 
the interesting feature of displaying a singularity in the 
spirit of multiscale expansion Q. As a result, meander- 
wavelengthes larger than the ones predicted by Bales and 



Zangwill can be explained, however, without reaching the 
order of experimental observations. Other efforts cen- 
ter around the investigation of an extra diffusion current 
along step edges . Such a current triggers an asymme- 
try in energy barriers for atoms attaching at kinks from 
different step edge directions, the so called Kink Ehrlich- 
Schwoebel effect (KESE) §. KESE can either destabi- 
lize or stabilize steps depending on whether the slope of 
the meander-instability is greater or less than one. In |^ 
calculations including a stabilizing KESE-current along 
steps were compared to the further thesis that islands 
might nucleate at step edges and thereby constitute a 
still different wavelength of instability. It appeared that 
the latter thesis fits best to the experimental wavelength- 
values. 

New light was shed on these investigations by exper- 
iments which revealed that a morphology qualitatively 
different from the meandering one can develop on vic- 
inal surface as well, with growth conditions being ex- 
actly the same as for the meandering-morphology ex- 
cept a different angle of miscut . It can be character- 
ized by the absence of one global growth direction and 
thereby resembles the so called degenerated morphology 
|l0| . Thus the insufficiency of theoretical approaches to 
explain the difference between the wavelength of instabil- 
ity predicted by Bales and Zangwill and the one observed 
experimtentally was extended by their failure to explain 
this degenerated morphology. Note that the surface of 
the in-phase-meandering morphology displayed in the 
STM topograph in Fig. 1(b) |^ was wrongly identified 
as Cu(l,l,17) and is in fact a Cu(0,2,24) surface The 
true Cu(l,l,17) on the other hand is the one to display 
a degenerated morphology 

In this letter I will introduce a novel model which 
can explain the rise of two different morphologies as 
well as their basic wavelengthes including KESE currents 
plus entropic step-step repulsion in the classical Burton- 
Cabrera-Frank (BCF) model [|l2|. It is the competition 
between the destabilizing and the stabilizing forces com- 
bined in this extended BCF model (EBCF model) which 
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results into two different kinds of morphologies depend- 
ing on the difference in magnitude of orthogonal driving 
forces. A careful analysis of this interaction also reveals a 
new basic wavelength of instability which is in the order 
of the experimental observations. 

To get an impression of the basic ideas underlaying the 
EBCF let us start with the one sided BCF model, i.e. 
the model formulated by Burton, Cabrera and Frank in 
1951 for the case of complete supression of attachment of 
adatoms to decending steps. This one-sided BCF model 
constitutes a moving boundary problem of a diffusion- 
relaxation equation for the dynamics of the adatoms on 
the terraces and two boundary conditions for the con- 
servation of mass and the conservation of energy at the 
steps, respectively: 
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Here c denotes the areal adatom density, Cg^ the equi- 
librium concentration for a straight step and k the step 
curvature. Equations (2) and (3) are to be evaluated at 
the front of each step. For vicinal Cu surfaces in the 
temperature range of the experiments under discussion a 
current involving the diffusion coefficient Dm along the 
kinked steps is operative Q. Following the notation of 
Pierre-Louis et al. Q this diffusion along steps plus its 
anisotropy due to KESE can be included in the mov- 
ing boundary problem above extending Eq.(3) by —dxJ, 
with J — Jk + Je {Jk, Je as in Eq.(2)-(4) in the pa- 
per of Pierre-Louis et al). Note that imlike the authors 
of j|] I consider 7 to be the step stiffness instead of the 
line tension as the precise measure of the step's tendency 
to straighten and thus reduce its kink density. As a 
consequence it is taken into account that just as sur- 
face undergoing a faceting instability develops into the 
surface-tension-minimizing morphology, instabilities at 
step edges have ideal orientations which can minimize 
line tension as well. This is crucial for the model since 
this is exactly the point in which the Cu(l,l,17) and the 
Cu(0,2,24) vicinal surface differ. Both surfaces are made 
up of (001) terraces of monoatomic hight and mean ter- 
race width / = 21. 7A. However, for the Cu(l,l,17) sur- 
face straight step edges before the growth of instabilities 
are parallel to the <100> direction, whereas they are 
parallel to <110> in the case of Cu(0,2,24). As a con- 
sequence the direction of minimal step stiffness, which 
can be determined to be in <130> via an embedded 
atom method (in analogy to [|l3| for surfaces), encloses 
an angle of appr. 63.5° with the straight steps in case 
of Cu(0,2,24), whereas the angle is appr. 71.5*^ in the 
case of Cu(l,l,17). Since growth proceeds in the direc- 
tion of minimal step stiffness, these are the orientations 



for the initial formation of the instabilities. They trigger 
destabilizing Jk currents of different magnitudes along 
the step edge forming in this direction. For to evalu- 
ate Jk a quantitative value for the KESE-length Lg (ref. 

Eq.(3)) is required. The relevant kink energy bar- 
rier is the one directed in <110>, giving 0.518eV us- 
ing an effective medium theory potential ||l^,|2^. The 
diffusion barrier for the jump of a single adatom along 
the step edge in this direction is 0.399eV |14[2^, lead- 



ing to a KESE-length Ls of 201.2. The resulting jk 
current densities for Cu(0,2,24) versus Cu(l,l,17) dif- 
fer by one order of magnitude: for Cu(0,2,24) jk takes 
a value of 7.118 • IQ-'^iA^s)-^, whereas at Cu(l,l,17) 
jk = 6.396 • 10"*(A2s)~^ Despite their different mag- 
nitudes in both cases KESE currents are destabilizing, 
favoring unsaturating amplitude growth. Wavelengthes 
of instabilities turn out to be even less then in the BCF 
model. As a consequence taking into account KESE as 
only additional driving force in the BCF model is not well 
suited to explain the experimental work by T. Maroutian 
and coauthors. An obvious antagonist of the destabiliz- 
ing KESE is the repulsion due to the succeeding step. 
Entropic just as well as elastic interactions ||l^ have to 
be taken into account. Since the step interaction en- 
gergy A considering an elastic or a dipole momentum 
is small compared to the entropic repulsion in the tem- 
perature ranges under discussion ||T^, it is sufficient to 
extend Eq.(2) by taking into account the surpression of 
step wandering: 
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The additional term {iTkBTYl&l 7 is the step interaction 
parameter with I equal the width of the terraces. Equa- 
tions (1), (2) and (3) constitute the extended BCF model 
(EBCF). Its simulation with T=280K, F=3-10-3 ML/s, 
1=21.7 A, Dm = lO-^cmVs, Ce, = 8.208 • lO^^A^i and 
7= 1.034-eV/A yields the following two morphologies for 
the Cu(0,2,24) surface (Fig. 1) and the Cu(l,l,17) (Fig. 
2), respectively. 

The four components regulating growth of instabilities 
in the EBCF are 

1. driving of growth via the gradient of the adatom 
diffusion field normal to the interface setting the 
length scale La for a primary wavelength of insta- 
bility 

2. restore via step stiffness (corresponding length 
scale: F) 

3. driving of amplitude growth via KESE (length 
scale: L^, ref. |@]) 
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4. restore via entropic repulsion (length scale: Ls = 




800 1450 2100 GU 



FIG. 1. In-phase meandering according to numerical sim- 
ulations of EBCF. The parameters of the simulation are cal- 
ibrated to Fig. 1(b) in GU referes to grid units of the 
underlaying numerical grid. A train of 15 steps with periodic 
boundary conditions in lateral direction is displayed. Hori- 
zontal boundary conditions are periodic as well. Space cali- 
bration leads to roughly 10 grid units corresponding to lA. 
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FIG. 2. The degenerated morphology as obtained in sim- 
ulations. Same as in Fig.^ except a different scale of IGU 
corresponding lA and taking into account the different vici- 
nal (namely Cu(l,l,17) versus Cu(0,2,24)) via a different di- 
rection of minimal step stiffness and resulting different mag- 
nitude of KESE (refer to text). A train of 25 steps with first 
step growing towards infinity is displayed. Steps feel the in- 
fuence of overlapping fields ||l^. This simulation result is in 
good agreement with the degenerated morphology recorded 
for Cu(l,l,17) after deposition of about 20 ML at F=5-10"^ 
ML/s at surface temperature 285 K via STM topography H. 



In the situation depicted by Fig. 2 the magnitude of the 
KESE current derivative dxJ is neglectable small. The 
remaining forces due entropic and step stiffness effects act 
in directions perpendicular to each other. The result of 
this competition of prependicular forces of equal strength 
is the degenerated morphology as observed in other con- 
textes with an analogous competition of driving forces 

In contrary in Fig. 1 the KESE current does contribute 
to the evolution of surface morphology. The magnitude 
of the driving force due to this current is of an order 10 
smaller then the magnitude of the adatom diffusion field 
gradient in direction of minimal step stiffness. The pre- 
cise factor is evaluated via simulation and reads 11.129 
for the simulation run done for Fig. 1 averaged over 500 
timesteps. Due to the dominance of La growth can pro- 
ceed in the direction of minimal step stiffness. 

For to understand the basic wavelength of the instabil- 
ity the two new length scales Ls and Lg, set by entropic 
interaction as in Eq. (2) and by an anisotrop diffusion cur- 
rent along the step as in Eq.(3) respectively, have to be 
taken into account. Eq.(l),(2) and (3) constitutes a sys- 
tem with a type I bifurcation |l^. It's dispersion relation 
can be evaluated in analogy to a dilute binary alloy un- 
dergoing directional sohdification |l^. For F <C the 
resulting expression in leading order reads 

A„ = V2n{ ^^-l-^' f'' + 2V27r^. (4) 

It replaces Eq.(l) in Expr.(4) is sufficient in the most 
interesting temperature range 2.8 to 3.8 x 10~^K~^. It 
displays deviations from an Arrhenius-type behaviour, 
which increase with increasing temperature. Neverthe- 
less it is in good agreement with the experimental data 
(Fig. 3). Simulations of the EBCF model support this 
expression as well. 
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FIG. 3. The '+'-data points arc taken from |5[ and give 
the wavelength measurements recorded there. The full line is 
the solution of (4). Crosses with errorbars are data points ob- 
tained from simulation with the EBCF. There are limitations 
to a precise wavelength-measure in the simulations. Due to 
the horizontal periodic boundary conditions wavelengthes are 
always a divider of the system's width. Enlargening this sys- 
tem width systematically until the rise of a further cell as 
well as restricting the width until the extinction of one of the 
cells allows to find upper and lower bounds. Data points are 
the mean of these bounds which themselves are taken into 
account via the error bars. Thus each data point is a result 
of up to 15 simulations with different system widthes. 

The extended BCF explains two relevant new results 
for growth at vicinal surfaces, namely the experimentally 
observed deviation from the Bales-Zangwill instability 
and the rise of a degenerated morphology. These re- 
sults are related to experimental findings in the group of 
H.-J. Ernst. The new basic wavelength predicted theo- 
retically from the EBCF-model can be observed as ini- 
tial perturbation of the step at Cu(0,2,24) just as well 
as Cu(l,l,17) surfaces. The qualitative difference be- 
tween their morphologies, i.e. whether a surface displays 
an in-phase-meandering (Cu(0,2,24)) or a degenerated 
structure (Cu(l,l,17)) after deposition of a few mono- 
layers, depends on the magnitude of the angle between 
the direction of destabilizing and restoring force. Time 
scales for the development from initial perturbation to 
the full degenerated morphology could not be compared 
with experiments due to the lack of experimental data. It 
seems an interesting question whether in situ transitions 
from one morphology to the other can be obtained via 
a change in the ratio of driving forces (e.g. by lowering 
temperature during the experiment). If this morphology 
transition displays features which resemble a true phase 
transition remains to be investigated. 

The diffusion-relaxation Eq.(l) is solved on a 
quadratic grid. The interface is discretized separatly by 
curvilinear segments and respective interpolations from 
the interface to the grid. Details of the scheme and its 
parallelization will be given in a further paper. 

I thank T. Maroutian and M. Rusanen for sending me a 
STM-picture of the true Cu (1,1,17) surface-morpology 
as well as a respective preprint prior to publication. Com- 
ments by H. Miiller-Krumbhaar to focus my attention on 
a precise understanding of the basic wavelength of insta- 
bility arising in the EBCF-model are gratefully acknowl- 
edged. Part of the computations were done on the T3E 
and the Origin 3800 of the URZ at Dresden Technical 
University. 

Note-The two surfaces{Cu(l,l,17) and Cu(0,2,24)) have 
also been modeled in on the basis of a kinetic Monte 
Carlo algorithm without taking into account entropic re- 
pulsion. The agreement with experimental findings re- 
mains restricted, nevertheless it is a very precise study 



of the different effects of KESE at the different surfaces. 
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